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Abstract 

We investigate the spin-independent part of the virtual Compton scattering 
(VCS) amplitude off the nucleon within the framework of chiral perturbation 
theory. We perform a consistent calculation to third order in external mo- 
menta according to Weinberg's power counting. With this calculation we can 
determine the second- and fourth-order structure-dependent coefficients of 
the general low-energy expansion of the spin-averaged VCS amplitude based 
on gauge invariance, crossing symmetry and the discrete symmetries. We 
discuss the kinematical regime to which our calculation can be applied and 
compare our expansion with the multipole expansion by Guichon, Liu and 
Thomas. We establish the connection of our calculation with the generalized 
polarizabilities of the nucleon where it is possible. 
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I. INTRODUCTION 



Recently, there has been greatly increased activity in the field of virtual Compton scat- 
tering (VCS) off the nucleon both on the experimental |I]-§J as well as on the theoretical 
f6HTT| side. At the c.w. electron accelerator MAMI (Mainz) the Al collaboration has taken 
initial data 0, and proposals for similar experiments have been developed at CEBAF (New- 
port News) fm or are being prepared at MIT-Bates ||. Compared with ordinary (real) 
Compton scattering, VCS off the nucleon offers a much greater variety of experimental pos- 
sibilities, because the response of the hadronic system to the electromagnetic probe can be 
investigated by independently varying the energy and the momentum of the initial state 
photon and thus allows for probing nucleon structure in both transverse and longitudinal 
modes over a wide kinematic range, covering both the perturbative and nonperturbative 
regimes. In the former one can describe the VCS process in terms of a perturbation 
series using the interaction vertices of QCD with quarks and gluons as explicit degrees of 
freedom. Studies of this perturbative regime of QCD require high resolution. Using an 
electromagnetic probe, this means that one has to perform an experiment wherein the mo- 
mentum transfer to the hadronic system mediated by a virtual photon is large, so that we 
are in the kinematic range of deep inelastic electron scattering. On the other hand, a long 
history of experimental facts suggests that the effective degrees of freedom of a hadronic 
system at low energies are not quarks and gluons but baryons and mesons or - in our special 
case - nucleons and pions. If one performs an electron scattering experiment which does not 
involve large momentum transfer to the hadronic system, one probes the confined phase of 
strongly interacting matter and does not resolve its underlying quark structure — our paper 
will deal with this nonperturbative phase of a hadronic system. (Another interesting line 
of work is probe the regime between the perturbative and nonperturbative phases. It is not 
clear at which momentum transfer such a "phase transition" will take place, and this will 
no doubt be addressed by future experiments. ) 

Despite much effort, even in the nonperturbative regime there remain many questions to 
be investigated in order to complete our picture of the nucleon and other hadrons. Important 
quantities which probe the compositeness of a system are its electromagnetic polarizabilities, 
which characterize the response of the system to an external electric or magnetic field |12| . 



In real Compton scattering, the electric and magnetic polarizabilities of the nucleon — ao 
and Pq — appear as the first model-dependent coefficients beyond the low-energy theorem 



(LET) of Low JT3[ and Gell-Mann and Goldberger [fL4|| , and have been measured both for 
the proton and the neutron JOJ. These polarizabilities have also been calculated within 
various models (for an overview see, e.g., |1|[I7|]). The LET of real Compton scattering 
has recently been extended to include virtual photons in P,|S[], and the generalizations of 
the polarizabilities were defined in ||, where also a first prediction for these "generalized 
polarizabilities" was obtained within the framework of a non-relativistic quark model. These 
results were refined in || to include recoil corrections. In [|7| an effective Lagrangian approach 
including nucleon and two-pion resonance mechanisms was used to determine the generalized 
electric and magnetic polarizability as a function of the initial photon momentum — a (| q |) 
and P (| q\). A very recent field theoretical calculation |L0| investigates the spin-independent 
polarizabilities within the framework of the linear sigma model. We will herein utilize the 
technique of heavy baryon chiral perturbation theory and compare our results with these 
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earlier calculations. 

It is well-known that the nonperturbative region can in general be successfully described 
in terms of the interaction of baryon and pseudoscalar meson degrees of freedom. In the last 
decade a new approach, chiral perturbation theory 0,0, was developed which describes 



the low-energy regime of QCD in terms of these effective degrees of freedom (pions and 
nucleons) while simultaneously requiring the symmetries of the underlying gauge theory 
and has yielded remarkable results. It was first applied to the sector of pseudoscalar mesons 
(see, e.g., |2(J for a pedagogical introduction) and then extended to case of pion-nucleon 



interactions The most recent version, which is known as the "heavy baryon formulation" 
of chiral perturbation theory [p^-2^|, uses techniques, which are well-known from heavy 



quark calculations. It allows a consistent power counting scheme to be developed, which 
was not possible with the former relativistic formulation. Real Compton scattering and 
the polarizabilities of the nucleon were among the first calculations to be performed in the 
nucleon sector of chiral perturbation theory [^,^-^], and the experimental verification of 
those predictions provides an important test of its validity. 

In this paper we extend these predictions to the VCS case and point out the strengths 
as well as the limitations of our approach. After a discussion in section II of the kinematics 
and the amplitude structure for VCS, we will in section III briefly touch upon the basic in- 
gredients of heavy baryon chiral perturbation theory and then proceed to calculate the tree 
and loop diagrams to third order in the chiral expansion with respect to external momenta 
- 0(p 3 ). In section IV we analyze the low-energy structure of the VCS amplitude, which 
follows from general principles like gauge invariance, crossing symmetry and the discrete 
symmetries pJTjJ, and predict the structure-dependent constants for the spin-independen^ 



part of the VCS amplitude, which are beyond the predictive power of such a general ap- 
proach. Finally, in a concluding section V we discuss the connection of those structure 
constants to the generalized polarizabilities defined in and relate the two approaches. 

II. KINEMATICS AND CHIRAL EXPANSION 

We begin our discussion of VCS off the nucleon by specifying our notation for the process 

7*(e", <f) + N(rf) - 7(^T q") + # (pjf) • (1) 

Here the nucleon four- momenta in the initial and final state are denoted by pf = (Ei,pi) 
and p^ = (Ef,pf), respectively. Since we will not in this paper discuss the component of the 
VCS amplitude which depends on nucleon spin, we have omitted indices for the nucleon spin 
states. The initial (final) state photon is characterized by its four-momentum = (u, q) 
(g'M = (u)',q')) and polarization vector = (s°,e) (e'^ = (e'°,e')). Whereas the final state 
photon is assumed to be real, 

q' 2 = ^ - q' 2 = , (2) 



The analogous spin-dependent calculation will be described in a subsequent publication. 
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the initial state photon is taken to be space-like, i.e. 



q i = -Qi<0. (3) 

Since our discussion refers to an electron scattering experiment, wherein the virtual photon 
is the exchanged particle between the electron and the hadronic current, we can write the 
polarization vector of the virtual photon as 

= eUe^^Ue/q 2 , (4) 

where u e and u e > are the initial and final state electron Dirac spinors, 7 M is a Dirac matrix 



(see, e.g., [J28| ) and e =| e |~ y An/ 137 > 0. Finally, we define the Lorentz invariant 
momentum transfer 

t = {q- q'f , (5) 

which will be useful in the following discussion. We evaluate the VCS amplitude in the cm. 

system, 

Pi = -<t, 

Pf = -q' , (6) 

and for our calculation we select a special frame of reference, 

g=(0,0,l), 

q = (sin 0,0, cos fl), (7) 

wherein we define an orthonormal set of basis vectors: 

e z = q, (8a) 
e y = q x q I sin 6 , (8b) 
e x = e y x e z . (8c) 

The complete VCS amplitude can in general be expressed in terms of three independent 
kinematical quantities, e.g., u' , \ q | and 9 ||. Note that electron scattering kinematics 
implies u' <| q \ which follows from energy conservation and Q 2 > 0. It turns out that our 
choice of variables facilitates the chiral expansion of the amplitudes, allowing chiral power 



counting ||18[ to be applied straightforwardly. The point here is that the chiral expansion is 
an expansion in terms of external momenta p, where in VCS on the nucleon the scale is set 
by the nucleon mass M and 4ttF. The pion decay constant has the value F = 92.4 MeV. 
Now it is a crucial point to define what is meant by the term "external momentum" . In the 
following analysis we will take the view that each component of the four-momenta of the 
photons and of the three-momenta of the nucleons (see Eq. (^)) has to be "small" compared 
to M.0 As a consequence we will count the terms u'/M and | q \ /M to be of the same chiral 



2 As a practical matter, since delta degrees of freedom are not included directly, the region of 
applicability is more appropriately to' , \q\ < m^. 
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order, whereas a term like | q \ 2 /M 2 is suppressed by one chiral order with respect to the two 
preceding terms. This has important consequences: Let us, e.g., consider the difference of 
the initial and final state photon energies, which occurs frequently in the subsequent heavy 
baryon calculation. From energy conservation one obtains 

where r stands for either u' or | q \ . From Eq. (||) we infer that we must count the difference 
of the initial and final state photon energies as being one order higher than both energies 
themselves. Now imagine that we have evaluated a diagram which is of the chiral order 
O (p n ) and that part of this amplitude is proportional to w - uj'. We find that this compo- 
nent of the amplitude can then be neglected because it contributes at the same order as a 
genuine O (p n+1 ) diagram. Specifically we will find below that the loop diagrams we wish 
to calculate are already O (p 3 ) and contain many pieces which are proportional to uj — uj' . 
Such contributions vanish from the beginning in the case of real Compton scattering, and 
we can also neglect such terms in our 0{p 3 ) VCS calculation, because they would contribute 
at the same order as a O (p 4 ) contribution, which is not addressed in this paper. Finally, it 
is useful to present the expansion of q 2 and t in terms of the three independent kinematical 
quantities, 

g 2 = ^ 2 -|g| 2 +0^ , (10a) 
t = -co' 2 - | q\ 2 +2oj' | q\ cosfl + O , (10b) 

which we will also use below. 



III. CALCULATION IN HEAVY BARYON CHIRAL PERTURBATION THEORY 

In this section we will extend the O (p 3 ) chiral heavy baryon calculation for real Compton 
scattering [P5,p? to the case where q 2 < 0. In the case of real Compton scattering in the 



forward direction, a O (p A ) calculation has been performed in [25|,26], which yields only small 



corrections to the O (p 3 ) result. Guided by this observation, we expect our O (p 3 ) calculation 
to provide a reasonable estimate in the kinematical region we are considering. The invariant 
amplitude for VCS can be written in the form 



Mvcs = -ie 2 E^M» = -*e V;M^ . (11) 

We will utilize the Coulomb gauge, 

e'" = (0, e") > e'-q' = (12) 

for the real photon. We decompose the space components of the virtual photon polarization 
vector into a purely transverse and a purely longitudinal part in our frame of reference, 
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e = e T + e L , (13) 

e T = e x e x + e y e y , (14) 

e L = e z e z . (15) 

Technically we calculate Aiycs with an initial photon polarization vector 

q 2 g 2 

a M =(0, a), a = e T + — e • gg = e T + — e 2 g (16) 



with g-a ^ | 29| , which simplifies the calculation by substantially reducing the number 



of Feynman diagrams which must be considered. The invariant amplitude can also be 
decomposed into a transverse and a longitudinal part. Using current conservation, 

q^ = 0, (17) 

Eq. ( |Tl~D can be written as 

M vcs = ie 2 (e T ■ M T + ^e z M z ^j . (18) 

The transverse part of the invariant amplitude consists of eight independent structures. We 
will use a notation similar to that given in [^7] , 

e T . M T = e'* ■ e T Ai + e'* ■ qe T ■ q'A 2 

+ io ■ (e '* x st) A 3 + ia ■ (q' x q) e'* ■ EtA^ 

+ io ■ (t* x q) e T ■ q'A 5 + ia ■ (e'* x q') e T ■ q'A 6 

- ia ■ (e T x q') t* ■ qA 7 - ia ■ (e T x q) t* ■ qA 8 , (19) 

where <7j {i G {x,y,z}) are the Pauli spin matrices. (Note that in the special case of 
real Compton scattering time reversal invariance imposes two additional constraints on the 
amplitudes, as a result of which we find A§ = A 7 and Aq = A$ .) Only two of these 
amplitudes — A\ and A 2 — are independent of nucleon spin. For the longitudinal component 
one finds four independent structures, 

M z = f ■ qA 9 + ia ■ {q' x q) t* ■ qA w 

+ia ■ t* x qA n + id ■ e'* x q'A 12 , (20) 

and in this case a single amplitude — A 9 — is spin- independent. In this paper, we restrict 
our consideration to this spin- independent component of the matrix element, 

M n v °cs spm = Mvcs - Mfos 

(21) 



ie 2 



e'* ■ e T A x + t* ■ qe T ■ q'A 2 + t* ■ q^e z A 9 



which can also be obtained from the full amplitude Mvcs as 



M non-s pl n = 1 tj . {Mvcs) (22) 
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We calculate M.™Q S spm using the standard chiral perturbation theory Lagrangian in the 
heavy baryon formulation to O (p 3 ) in the nucleon sector p3|j2^| , 



C 



ttN 



(i) 

ttN 



r m , as) 



N i 



with 



C 



(i) 



nN 

(2) 
wN 



£ (2) 



N v (iv-D + g A S-u)N v , 

(v ■ D? 



1 -NAD-D 



2M 



2 c [iv pa v u 



/r(l + 4c 6 ) + 2^ s ^(l + 2c 7 



As) 



2M 2 



X £ 



fiupcr 



S a iD p + h.c. AT, 



(23) 



(24a) 



(24b) 



(24c) 



where £0123 = 1- We have only kept those terms which contribute to a O (p 3 ) VCS calcula- 
tion. In particular terms linear in the photon fields, which vanish in our gauge, have been 
omitted. Moreover, we note that contributes only to the spin-dependent piece of the 
VCS amplitude and is irrelevant for the spin-independent part. The velocity-dependent field 
N v is projected out from the nucleon Dirac spinor ^n. 



N v = exp [iMv ■ x] P+*jv , 
where the projection operator is given by 

where v M is a velocity vector satisfying v 2 — 1. The covariant derivative is defined as 



D, = d, + T,- i„ W 



where 



- {u ] (dp - irj u + u(d^- iZ M ) m 1 "} , 



" 2 
u = , 

[/ = exp (ir • 7r/F) 



(25) 



(26) 



(27) 



(28a) 

(28b) 
(28c) 



In Eq. Q28c|) r are the conventional Pauli isospin matrices, while 7? represents the interpo- 
lating pion field. The field strength tensors are defined as 



,(«) 



fJ,U 



z[rV1 



jjf = 0T - d u F - % [F, r 



(29a) 
(29b) 
(29c) 
(29d) 
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and to the order considered the corresponding coefficients are related to the anomalous 
magnetic moments, 



c 6 
c 7 



Hv 1 _ 1 / 
4 4 ~ 4 [fip 



2 



(flp + fir, 



(30a) 
(30b) 



In our application the right- and left-handed currents, and Z M , are the isovector part of 
the electromagnetic current and vff is the isoscalar piece, 



6 A 

e . 



6 -A 



The Pauli-Lubanski spin vector is given as 



(31) 



(32) 



with o^ v = | [7^, 7"]. Finally, we define the quantity 

= % {V {dp - ir M ) u u (d^, Up) w f } . (33) 
In the heavy baryon calculation we can write the initial and final nucleon momenta as 



Pi 



where if and tj denote off-shell momenta. We choose v M 
where the polarization vector a M is given in Eq. (pf). 

In the pion sector we can restrict ourselves to the O (p 2 ) Lagrangian, 



(34) 

'1, 0, 0, 0) such that v ■ a = 0, 



r(2) 



-tr 



where 



(35) 



(36) 



represents the covariant derivative. Here we have omitted the usual mass term, because 
we only need to consider pion-photon vertices generated by C^}. We also observe that 
power counting arguments, in principle, require the inclusion of a tree diagram with a tt° 
in the i-channel in our calculation (Fig. [I] (d)). For this diagram we need the lowest order 
Wess-Zumino Lagrangian p0|j31[| , which has odd intrinsic parity and is O (p 4 ). However, this 
diagram does not contribute to the spin-independent part of the VCS amplitude, A4^g Sptn , 
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but only to its spin-dependent counterpart, Mycs- For this reason we will not discuss it in 
the following. Finally, we wish to emphasize that we do not require any additional diagrams 
compared to the calculation for real Compton scattering |23| . The complete set of diagrams 
we have to calculate is given in Figs. [I] ((a) s-channel, (b) w-channel, (c) contact diagram) 
and ^ (loop diagrams). In the following we will treat the tree and loop parts of the amplitudes 
separately, 

Ai = A l r + A[ oop . (37) 



A. Tree Diagrams 

We expect the s- and it-channel Born diagrams (Fig. [I] (a,b)) and the contact diagram 
(Fig. [L] (c)) to generate the structure- independent part of the amplitude which was predicted 



for real Compton scattering by Low [y| and Gell-Mann and Goldberger [14]. In h has 
been shown that the heavy baryon approach reproduces this LET, which must be the case 
for any consistent calculation based on gauge invariance, Lorentz invariance and crossing 
symmetry. For the case of VCS we obtain 

AT = ~ (1 + n) , (38a) 
K" = l^f 2 ^+r 3 )\i\, (38b) 

' 4r " = -2S (1 + T3) + p (1 + T3)l « |c08 '' 

Recently the LET has been extended to the case of virtual Compton scattering on the nucleon 
HH]. We find that our results at O (p 3 ) exactly agree with the predictions, which behave as 
1/M and 1/M 2 in ||. However, the terms of order 1/M 3 which have also been derived in || 
as well as the q 2 corrections to the zeroeth order nucleon form factors, which are proportional 
to (4kF) 2 > P^l are beyond the predictive power of our O (p 3 ) chiral calculation. [] We note 



that we do not generate any contributions to the spin-independent part of the VCS amplitude 
in the case of VCS on the neutron, because the photon only couples to the nucleon charge 
at tree- level in the spin-independent sector, but not to its magnetic moment. 

B. Loop Diagrams 

In contrast to the three diagrams in Figs. [L] (a) - (c) discussed above we expect the 
loop diagrams in Fig. |] (and the t-channel diagram from Fig. (d)) to contribute to the 



'Such effects will be considered in a subsequent 0(p 4 ) publication. 



9 



"structure-dependent," and thus model-dependent terms beyond the LET. Qualitatively, 



this can be understood as follows: In the sense of Low's method 32 the s- and w-channel 



diagrams of Figs. [I] (a) and (b) generate the most singular terms, whereas the contact 
interaction of Fig. [I] (c) is required by gauge invariance. On the other hand, the one-particle 
irreducible loop diagrams of Fig. ^ yield regular contributions to the VCS amplitude, which 
are of higher order in the external momenta. These contributions are not predicted by gauge 
invariance alone. Nevertheless, we will see that an extended low-energy expansion will enable 
us to derive some constraints for a calculation of the structure-dependent amplitude (see 
section [TVp . 

In the calculation of the loop contribution to the amplitudes, A l ° op , it turns out, as 
already pointed out above, that the number of loop diagrams is reduced to 9 by the choice 
of the gauge and of the velocity v. By power counting arguments, the only possibility 
of generating O (p 3 ) loop diagrams consists in using interactions from c£' N and C^J. As a 
consequence the photon- nucleon interaction vanishes at lowest order, C^ N = 0, and we need 
not consider any loop diagrams which contain this kind of interaction. Using the interactions 
of Eqs. (|24a|) , (|24b|) and (B^), we obtain the invariant amplitudes for the 9 diagrams in Fig. 



0, M.ycs t° M-vcs- The exact results are listed in the appendix . For some diagrams we 
cannot carry out the integrations over one or even two Feynman parameters analytically and 
could, in principle, proceed to evaluate these integrals numerically. However, since we want 
to establish a connection with the general low-energy expansion of the structure-dependent 
part of the VCS amplitude, we expand the expressions in Eq. ( |X2| ) in terms of the two 
external momenta, oj' and | q \, using Eqs. (|D|), ( |10a| ) and ( |10b| ). This expansion will scale 
with the mass of the particle propagating in the loop. Consequently, we will obtain a power 
series in r/m n where r stands for oj' or | q |. 4 We could in principle extend the r/m n 
expansion to an arbitrary given order within the framework of a O (p 3 ) calculation. In this 
context we want to point out that it is crucial not to confuse the expansion in r/m^ with 
the chiral expansion in r/M. The result of the chiral O (p 3 ) calculation, expanded to the 
order r 4 then reads 



j^loop 



F 2 nm 
4 



/2 / I -* I n 

— oj H oj q cosu 

96 192 1 H 1 



w + tt^—u | g | cos 6> 



1920 ml 
1 1 



384 ml 
1 1 
+ 960m2 

F 2 nm^ 



12 I - 1 2 

oj \q\ 



1920 ml 
1 1 



320 ml 



OJ 



q | 2 cos 2 



oj' | q \ 3 cos 6 



oj | q 



192 



384 m: 



(39a) 



-oj' 3 I q 



4 This expansion does not increase the chiral order, because r and are both treated as "small" 
parameters. 
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320 ml 
9l 1 



^ I H 

5 



960 ml 



F 2 7rm Tl 

7 1 

+ — - oo 



-OJ 



12 



17 1 



96 

'" 1 q | cos# 



-OJ 



1920 ml 

7 1 



/2 I -\2 

oj \ q\ 



960 ml ' 960 
In order to obtain v4g 00p , we made use of Eqs. fl8a|) and (|16). 



(39c) 



IV. LOW-ENERGY EXPANSION, STRUCTURE COEFFICIENTS AND 

POLARIZABILITIES 



In the previous section we saw that the O (p 3 ) heavy baryon result reproduces the terms 
required by Low's method in the tree-level amplitudes. This section will deal with a general 
parametrization of the structure- or model-dependent part of j\4™^g Spm . In JTTj] a general 
low-energy parametrization O of the structure-dependent amplitude for virtual Compton 
scattering on a spin-zero target, e.g. a pion, has been worked out. The corresponding 
expression for the nucleon is more complicated, as in Dirac space it involves a general 4x4- 



matrix. However, we can apply the parametrization of [|TT|] to the spin-independent part of 
the VCS amplitude off the nucleon. Indeed the form of the scattering amplitude for the 
spin-independent part of the reaction ([]]) is the same as for the process |33| 
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j(e»*,q'») 



+ ir(p 



(40) 



Before discussing the structure-dependent terms, one has to specify which convention has 
been used for splitting the total amplitude M.ycs into a Born contribution and a structure- 
dependent part. This issue is addressed in P,|TTf in quite some detail. We follow the 
convention of fHH where the Born contribution is calculated using Dirac and Pauli form 
factors, F\ and Fi- In fact, to the order considered here, the tree-level diagrams of Fig. 
HI (a) to (c) generate the same contributions to A\, A 2 and A$ as a covariant Born term 
calculation involving Fx and F 2 ||. 

The most general spin-independent contribution to the structure-dependent VCS ampli- 
tude can then be written as 



O = e^O^e'* = (1) + (2) + (3) + (4) + 0{k 5 ), 
where the terms of increasing orders of k (k stands for q or q') are the following: 
O« = 0, 

O^=g [e-q'e'*-q-q-q'e-e'*] 

+ci [(q + q ) ■ (pi + p f ) (e ■ (p, + p f ) e'* ■ q + q' ■ e (p, + p f ) ■ e'*) 
-2q ■ q'e ■ (pi + p f ) e'* ■ (pi + p f ) - 2q ■ ( Pi + p f ) q' ■ ( Pi +p f )e- e'*] , 

O^=0, 

(4) = g 2a q ■ q + g2bq 2 + 4g 2c q ■ (Pi + Pf) q' ■ {Pi + Pf) [e ■ q'e'* ■ q - q ■ q'e ■ e' 



(41) 



(42a) 



(42b) 
(42c) 
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+ [c 3a g • q' + c 3b q 2 + Ac 3c q ■ (p { + p f ) q' ■ (p { + p f ) 
x [(q + q') ■ ( Pi + p f ) (e • ( Pi + p f )e'*-q + e- q'e'* ■ ( Pi + pf)) 
-2q ■ q'e ■ (pi + p f ) e'* ■ fp { + p f ) - 2q ■ (p { + p f ) q ■ (p { + p f ) e ■ e'*} 
+c 3 [2q 2 (q ■ (pi + pf) q ■ (p { + p f ) £ ■ e'* - q ■ q'e ■ (pi + p f ) e'* ■ (p { + p f )) 

+ (q + q) ■ (pi + Pf) q 2 (e • (pi + Pf) s'* ■ q - e ■ q'e'* ■ (pi + p f )) 
+2 (q + q') ■ ( Pi +p f )q- q'e ■ qe'* ■ (p { + p f ) 
-Aq ■ (pi + p f ) q ■ (pi +p f )e- qe'* ■ q] . 

Here we have used the same notation for the unknown structure coefficients as in ||11|| . 
We then match our chiral calculation of the O (p 3 ) spin-independent amplitude with the 
low-energy expansion, Eq. (|4l|), demanding 



(42d) 



-ie 2 



i Aoop,non— spin 
JVl VCS 



(43) 



A consistent matching procedure requires that we make use of identical approximations for 
the kinematical quantities (see Eqs. (Q), (|10a|) and (|10b| )) in both the chiral loop calculation 
and the low-energy expansion, Eq. (f4lP . To be specific, we have neglected corrections of the 
type u'/M and | q \ /M (and higher). After this procedure we are able to determine the 
unknown structure coefficients in Eqs. ( |42b|) and ( |42dj) from our chiral calculation. Using 
Eq. ( pTTP we obtain the following parametrization of the spin-independent amplitudes: 



structure 



-go 



8M ci + uj' | q | cos 9g 



+UJ 



-92a - 92b - 16M^ 2c + 8M' z c 3 - 8M% a - 8M z 5 3b - 128M 4 £ 3 , 



r2~ 



r4~ 



+oo' 3 | q | cos 6 2g 2a + 92b + 16M^ 2c + 8M' z c 3a 



2~ 



^structure 



a structure 
A 9 



+lu' 2 | q\ 2 [g 2b - 8M 2 c 3 + 8M% b \ + J 2 \ q \ 2 cos 2 6 [-g 2a 
W | q | 3 cos 9 [-g 2b ] , 
' ' ~" ' r -9o] 



u | q 
+u' 3 

12 



q 



-92a 



g 2b - 16M 2 g 2c +u' 2 \q\ 2 cos 9g 2a + J \ q 



92b 



(44a) 
(44b) 



to 



+UJ 



/l 



- 8M 2 d lj 

-92a - 92b - 16M 2 g 2c + 8M 2 c 3 - 8M 2 c 3a - 8M 2 c 3b - 128M% C 



+uj' 3 \q\cos6 g 2a + 8M' A c 3a 



2~ 



+ uj \q 



g 2b + 8M' z c 3 + 8M 2 c 3b 



(44c) 



Demanding the validity of the matching relation Eq. fl43|) we obtain a system of linear 
equations for the structure coefficients, yielding 



go 



g2a 



i g\ 



192 F 2 vrm 
11 1 



g\ 



192 8M 2 F 2 vrm, 
1 



i g\ 



320 F 2 irml 



(45a) 
(45b) 
(45c) 
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92b = — ~4^T (45d) 
fe 1920 16M 2 F 2 vrm3 ' 1 ' 

c = _^ i_g\J_ f45n 

3 1280 8M 2 F 2 vrm3 ' 1 } 

\ \ g\ \ 

° 3a ~ 240 81^^^4' (45g) 

36 256 8M 2 F 2 Trml ' 1 J 

9 1 g\ 1 

C3c ~ 640 128M 4 F 2 7rm3 ' ^ lj 

Moreover, we wish to stress that it is actually a prediction of the low-energy expansion, Eq. 
( |44a|) , that the expression for A\ cannot involve the structures | q | 2 and | q | 4 . Similarly 
the expansion of the amplitude A 9 , Eq. ( |44cj ), excludes structures of the type to' \ q \ cos 9, 
uj' 2 | q | 2 cos 2 #, uj' | q | 3 cos# and | <f | 4 . When we compare these predictions with the 
result of our loop calculation, we find that each of these constraints is satisfied. We also 
note that our O (p 3 ) calculation gives a zero result for structures with an odd power of 



r (r G {a/, | q |}), which is required by the general expression in |T| obtained by imposing 
gauge invariance, crossing symmetry and the discrete symmetries. 

It is useful at this point to interpret the structure coefficients which we have obtained. 
The coefficients which originate from — go and c\ — are closely related to the electric 
and magnetic polarizabilities of the nucleon in real Compton scattering, cto and (3q. This 
can be seen by considering the limit of real Compton scattering in the cm. system, i.e. 



q\=uj = to'. If we apply the definitions for the polarizabilities, (see, e.g., [f27]1 ) 



«o + A) = — -rr^At (uj = 0, 9 = 0) 



■2 ^2 

8nduj 2 

/3o = -£(^)(^ = 0,* = 0), (46) 




to Eqs. (|39a|) and ( |39b|) , we determine 



a 



5e 2 g 2 A 



384vr 2 m T F 2 



A, = n Jf A m , (47) 
7687r / m 7r i' 2 



which coincide with the results found in [p3fl . We can then express the structure constants 
g and C\ in terms of a and (3 , 

An 

9o = —Po , (48a) 
e z 

5i = -^(«o + A)). (48b) 
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An alternative way by which to derive this relation without explicitly taking the real Comp- 
ton limit is to compare our results for A\, A 2 and A 9 with the low-energy expansion of the 
VCS amplitude of ||. We find that this method yields identical results for the polarizabilities 
«o and (3q once the transformation to our convention is made. 



V. COMPARISON WITH OTHER CALCULATIONS 

We now investigate the connection of our results with the multipole expansion of ||, 
where possible. We begin by recalling the primary features of such an expansion. Guichon, 
Liu and Thomas parametrize the structure-dependent part of the VCS amplitude in terms 
of essentially the same kinematical quantities as do we, namely, uj', | q\, 9. However, they 
expand the structure-dependent part of the amplitude in terms of uj', keeping only the first 
non-vanishing, linear term. (They do not consider uj' 2 terms because they are of higher order 
in their expansion.) The multipole expansion then generates various combinations of cos# 
and powers of | q | and, furthermore, suggests the definition of generalized polarizabilities 
characterized by angular momentum quantum numbers of the respective partial waves. The 
crucial difference in comparison with the low-energy expansion introduced above consists in 
the feature that the starting point of the Guichon et al. analysis is the uj' expansion around 
uj' = 0, wherein | q | can be chosen arbitrarily, (i.e. not necessarily small), once uj' has a 
small value. One can see the difference between the two expansion schemes most easily by 
looking at the quantity uj. Rewriting Eq. (|]), we obtain 

On the other hand, Guichon, Liu and Thomas use 

«U^-^+0(JlC). (50) 

From Eqs. (|49|) and ( |50l) we conclude that we will not, in general, be able to express the 
results of our calculation for the structure coefficients in terms of the generalized polarizabil- 
ities. In our approximation scheme we take the quantities uj and uj' to be equal and consider 
terms up to uj' 4 whereas only keeps one power of uj' and starts the expansion of uj with 
a term which due to its 1/M suppression would be higher order in our scheme. Using Eq. 
fl50|) the expansions of the mo del- dependent parts of the spin-independent amplitudes in Q 
read 



A non-Born = _ - / | q | ^ ^(11,11)0 (| -|) 

V 8 

-f r 'uj\^ P^°(\q\) 

-\uj'\q\ 2 P^{\q\)+0(uj' 2 ) , (51a) 

A non-Born = Jlu'\$\ p(H.H)0 (| g |) + O (uj' 2 ) , (51b) 



14 




A non-Born = 1*^ ^ p(01,01)0 (| ?)) + Q ^ ? (51c) 



where we had to modify the definition of Ag on ~ Born in §, because therein the Lorentz gauge, 
£ ■ q = 0, is used, while we work in the gauge defined in Eq. (|16|). Note that the generalized 
polarizabilities p( 11 ' 11 )° ; p(oi,oi)o anc j p(oi,i)o are f unc tions of | q\ only. 

Let us now compare Eq. ( |39aj ) with Eq. ( |51a| ). First of all, we note that in virtual 
Compton scattering the term containing ao is proportional to ujuj' before we have made any 
kinematical approximations. This can be seen from Eqs. ( [42 b| ) and ( [48 b| ) . When we insert 
Eq. (|49|) for uj the leading term will be proportional to uj' 2 , whereas applying Eq. (|50|) the 
leading term will be proportional to uj'uj | w / = o. The important observation here is that in 
both expansion schemes these terms have the same coefficient which is proportional to a^. 
We conclude that our result for the uj' 2 term in Eq. fl39a|) serves as a prediction for a in || , 




P <w.oi)o ( |7| =0)= - /* 0= *f g Q + 8 M%) . (52) 



'2 

With an analogous chain of arguments we find that the term uj' \ q | cos 9 in Eq. ( |51a|) can 
be directly compared with the corresponding term in Eq. (|39a|) , yielding 

pCiiADo ( | 0) = -^y|/3 = -y|, . (53) 

We see that the result for the magnetic polarizability (3q obtained from Eqs. ( |39a| ) and (|51a]j 
is consistent with the fact that both Eq. ( 39b ) and Eq. ( 51bf) can also be parametrized 



in terms of this polarizability as a pre-factor of the uj' \ q | term. As regards Eq. ( |51c| ) 
we observe that the parametrization in terms of p( 01 > 01 )° is consistent with our result for 
the pre-factor of the uj' \ q | term in Eq. ( P9"c| ). Now what about P^ 1 . 1 ) (| q |= 0) in Eq. 
( |51a|) ? In order to determine this polarizability from our calculation, one would have to 
calculate the coefficient of the uj' \ q | 2 term in Eq. (|39a|) . Due to the approximations of 
Eqs. (0), ( |10aj ) and ( |10b| ), a O (p 3 ) heavy baryon calculation does not generate a term of this 
type. So we cannot determine this polarizability in our calculation. We expect the leading 
term contributing to pC 31 - 1 ) (| q |= 0) to be at least 1/M suppressed. Hence its evaluation 
requires a consistent O (p A ) calculation, which is in progress. This is borne out by the linear 
sigma model calculation ||TD|| which finds that the coefficient of uj' \ q | 2 is indeed 1/M 
suppressed. Furthermore, in [[Hj it was shown numerically that p(° 1 ' 1 ) is not independent 
of the generalized polarizabilities a and f3. 

Let us turn to the remaining seven structure constants, g2 a , g2b, 92c, C3, C3 Q , c^b, cs c , 



which parametrize the r 4 terms in Eq. (|44|) . First we observe in Eqs. (|39a|) and (|51a|) that, 



strictly speaking, we will only be able to compare terms which contain uj' (and under certain 
circumstances also uj' 2 ) but not the terms with uj' 3 and uj' 4 . These are higher-order corrections 
in Eq. (|51aD but not in Eq. (|39a|), where uj' and | q\ are counted as being of the same order. 



This suggests that the multipole expansion of Eq. (|5T|) and the chiral expansion of Eq. (39) 
have different domains of application. We will discuss the appropriate kinematical region 
for each of the two expansions later. We conclude that the terms uj' 3 ,uj ,4: belong to the 
model-dependent piece of the VCS amplitude, but do not fit into the expansion scheme in 



15 



terms of the generalized polarizabilities of Eq. (0). Let us now discuss the uj' 2 \ q\ 2 cos 2 9 
term in Eq. fl39a|) . We cannot expect to find an analogous term in Eq. (|51a|) , because the 
quadratic dependence on cos 9 indicates a term of higher multipolarity which has not been 
considered in the special application of Eq. (|51a|) , where the angular momentum of the final 



state photon is fixed at L' — 1. As regards the uj' 2 \ q \ 2 terms in Eqs. ( |39a| ) and ( |39c| ) we 



find that the factors uj' 2 originate from the product ujuj' before we apply Eq. (H). This can 



be seen from Eqs. ( |42d[ ), fl44ap and ( |44c|) . As a result we might be tempted to interpret 

\- \~\ r\ rus /-\ fTi ni nnf T~i /~\ i I I 2, -f- / \-v > -v y i n r\ n s~\, ( I /~t I (~\ \ T-T/-nir/""\'T r^v urn -fi n s~\ /~\ ^ -fl-V^-v /~\vi -f- nj-iln^-i /"in in 



the coefficient of the uj \ q\ terms as ^=p a (| q |= 0). However, we find different solutions 



from matching the transverse and the longitudinal results of Eqs. ( |39a|) and (|39c|) with Eqs. 

Inspecting Eqs. ( f44a| ) and ( f44c| ) we observe that the difference between 



( plaD and ( pel) 

the transverse and the longitudinal result can be explained by the fact that the structure 
constant C3 enters the terms with a different sign. If we recall that the transverse amplitude 
A° op contains contributions from higher multipoles than L' = 1, we come to the conclusion 
that the uj' 2 \ q \ 2 term can also get contributions from such multipoles, in particular from the 
constant part of the Legendre polynomial for L' = 2. For this reason we cannot determine 
^4pa(| q |= 0) from this amplitude. In the longitudinal amplitude A l g° p , however, we find 
no contributions from higher multipoles. For this reason we use the longitudinal amplitude 
for the determination of (| <f |= 0) and, using Eqs. (|39cj) , ( |44cj ) and (|51cj) , arrive at 



d 



d I q 



-a(\q\ = 0) = |- (g 2b + 8M 2 c 3 + 8M 



2 c 3 b 



38407r 2 m 7r F 2 



(54) 



The situation is more obvious for -Ji^ -/3 (| q \= 0). If we compare the term uj' \ q\ 3 cos 9 in 
Eq. ( |39al ) with Eqs. (pa]) and (|5Ta]) we find 



d e 2 

A 1 ->\ 2 fi (I <? 1= °) = -7-^26 = „„ ... > , • 
d I q Y 47T 38407r^m 7r i' 1 



(55) 



From Eq. ( |51b|) we can read off that we could expect the same coefficient for the uj \ q | 3 
term in Eqs. (|39b| ) and (|44b|) , and we find that our results are consistent with Eq. (|5 lb[) . 
This is supported by the fact that we have no indication of higher multipoles entering 
in our calculation which could obscure the V = 1 results. 

Summing up, we have established a connection between the multipole expansion M and 



the low-energy coefficients [|TT|] go, ci, g 2 b and the linear combination C3+C3b, which is reflected 
in the following expansion of the electric and magnetic polarizabilities of the nucleon, 



a (I q\) = a 1 




50 m 



Air 



(- [g + 8M 2 5 1 



+ I <f I 2 (92b + 8M 2 c 3 + 8M 2 5 3fe ) + O 




(56a) 
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P(\q\) =0o(l + ' 




5 m 



^U-\?\ 2 ^o(m). (56b) 



We wish to point out Eqs. ( |56a| ) and ( |56bQ exactly agree with the corresponding terms 



in the relativistic field-theoretical calculation within the linear sigma model |[10|| . In com 



parison with the results of [§,0,0 our value for the electric polarizability of the proton, 
«o = 12.8 x 10~ 4 fm 3 is much larger and our value for the magnetic polarizability of the 
proton, p = 1.3 x 10 _4 fm 3 is smaller than in those calculations. (We have used the numer- 
ical values = 135 MeV, F = 92.4 MeV and = 1.26.) The slope of the corresponding 
generalized electric polarizability, (| q\= 0), is found to be considerably larger than in 

the effective Lagrangian and the constituent quark model |6]|| calculations. The slope of 
the generalized magnetic polarizability, -J^P (I Q |= 0), even has a different sign compared 
with other calculations. For the neutron we find the same analytical expressions for these 
quantities as for the proton. 

We cannot give an interpretation of the structure coefficients g 2a , g2c, c 3a , c 3c and the 
missing linear combination C3 — in terms of generalized polarizabilities. However, all 
structure coefficients in Eq. (|42d|) , not only the subset which we could interpret in terms 



of generalized polarizabilities, have to be considered if one performs an experiment in a 
kinematical region where us' is comparable with | q\. Furthermore, we find that, in addition 
to the convergence radius of the chiral Lagrangian, which is given by r/M, the inelastic 
threshold of single pion production sets a primary limit to our calculation. The low-energy 
parametrization in terms of structure coefficients is only valid in a regime below the pion 
threshold, where u and | q | do not approach m n . This is clearly a different kinematical 
regime than the one of the multipole expansion ||, which is predominantly suited for ar- 
bitrary I q I, which are much larger than the very small u'. (Of course, in principle our 
calculation is exact for to' < m n and | q |< m v if we keep the finite integrals in Eq. ( |A2|) 
instead of expanding the integrands, as done above.) 

In conclusion, we have investigated the spin-averaged amplitude of virtual Compton 
scattering within the framework of heavy baryon chiral perturbation theory to O (p 3 ). We 
have determined the generalized spin-independent electromagnetic polarizabilities of the 
nucleon |J. In order to obtain analytic expressions, we have restricted our calculation to 
cto and (3q of real Compton scattering and the slopes of the generalized polarizabilities with 
respect to | q | 2 . We have also performed an alternative low-energy expansion of the VCS 
amplitude JTTJ] which is not restricted to first order in the energy of the outgoing photon as 
is the multipole expansion. We made a prediction, based on chiral symmetry, for the 9 a 
priori unknown structure coefficients characterizing the spin-independent part of the VCS 
amplitude up to 0(r 4 ). All predictions of the loop calculation are the same for the proton 
and the neutron at 0(p 3 ), but will differ in a 0(p 4 ) evaluation. 
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APPENDIX A: LOOP AMPLITUDES 



Using standard techniques we obtain the following invariant amplitudes for the Feynman 
diagrams in Fig. |2] with the Lagrangian interactions, Eqs. (|24a| ) , ( |24fe|) and (|35|): 

M% S = -||f / J^yi N ^* ■ \. { l + ti \ q)+l + P- m l + %0+ N ^ ( Al& ) 



2F 2 J (2tt) c 

x \ \ m (Ale) 

v(l + U + q-q')+iO+P-m 2 n + iO+(l-qi) 2 -ml + tO+ '" 



M&s = ^J-^- d N v e>*.la.(2l + q ) 



2F 2 J (2tt) c 

1 11 

x T-, n 7 To n ; 9 N v, (Aid) 

v(l + t i + q-q') + i0+ l 2 -m 2 n + i0+ (I + q) 2 - m 2 + t0+ 



M% 8 = §| / T^zAV* • la ■ (21 - g) 



2F 2 7 (2tt) 

x \ \ (Ale) 

y{l + t i )+iO+P-ml + iO+(l-q) 2 -ml + iO+ V 7 

= §| / ■ Is'* ■ (21 + q') 

1 1 1 

x — t, s n ; 5 AT W , (Alf) 

v{l + U)+iO+l 2 -ml + iO+ (l + q f -ml + t0+ V ; 

M vcs = e ^J -^r S N v (v.(l + q'-q)v-l-l-(l + q'-q)) 

xe'* ■ (2 (l-q)+q')a-(2l-q) 

x 1 1 1 1 N fAl 

X v ■ (I + t~) + 20+ I 2 - ml + i0+ (/ - qf -ml + i0+ (I + q' - g) 2 - m 2 + i0+ ( lg) 

M vcs = ^0J -^- d N v {v-{l + q'-q)v-l-l-{l + q'-q)) 

xa-(2(l + q')-q) e'* ■ (21 + q') 

x — — ~r — \ : tt^- ^ N v , (Alh) 

v(l + t i )+i0+l 2 -ml + t0+(l + qf-ml + i0+(l + q>-q) 2 -ml + i0+ V ' 

MP cs = -^ J ^ d N v (v ■ (l + q> - q)v -l - l- (l + q' - q)) 

xe'*-a — — \ : i N v . (Ali) 

t;-(/ + t i )+i0+/ 2 -m 2 +z0+(/ + g'-g) 2 -m 2 +z0+ 
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We have only retained the spin-independent parts of the ci-dimensional integrals. We note 
that we do not find any dependence on the nucleon isospin, which means that the loop 
contributions for proton and neutron are identical. Carrying out the integration over the 
loop momentum I we find that the sum of all amplitudes is finite. The individual amplitudes 
can be written as 



M {1) 



vcs 



M 



(2) 
VCS 



~2~F~ 2 ~ 

% e 2 g\ 



Jo (u, ml) , 



Mf cs = 
M% s = 



Mf cs 



M% s 



2 F 2 a - £ '* J o (-w'.mj) , 

2^2 „i 

e'* / dxJ' 2 (-u/ (l-x) + u, ml) , 

J 

1 



i£9a 

F 2 
ie 2 g\ 



F 2 



ie 2 g 2 A 
F 2 

1 F 2 



a ■ e'* / dxJ' 2 (—a/ + u (1 — x) , ml — q 2 x (1 — 

J 

+e'* -qa-q [ dx- (2x - 1) JL (-J + u (1 - x) , ml - q 2 x (1 - 
Jo 2 v 

J? 2 A a ■ e'* J dxJ' 2 (ujx, ml — q 2 x (1 — x)) 

+e'* -qa-q J dx^ (2x — 1) J' G (ujx, ml — q 2 x (1 — x)) 

a ■ e'* j dxJ' 2 (—uj'x, ml) , 

/ dx dy 
Jo Jo 

-5(1-1/) 



x 



a ■ e 



x J'l (u (y + x (1 - y)) - uj'y, ml - t (1 - x) y (1 - y) - q 2 x (1 - x) (1 - yf) 
' 't - q 2 ) (1 - yf (y + x (1 - y) + y (1 - x)) 



+ 



+q 2 (1 - x) (1 - yf (y + x (1 - y)) - J 2 y (1 - yf 

-lu 2 (1 -x)(l- yf (y + x(l- y)) + ujuj' (1 - yf (y + x (1 - y) + y (1 - a;)) ) 

x J'l (w(y + x(l- y)) - u'y, ml - t (1 - x) y (1 - y) - q 2 x (1 - x) (1 - yf) 

+a ■ q'e'* ■ q [ ((1 - x) (1 - y) 2 (1 - 6y) + y (1 - y) (x + y (1 - x))) 
xJl(u(y + x(l- y)) - tv'y, m 2 w - t (1 - x) y (1 - y) - q 2 x (1 - x) (1 - y) 2 ) 
+ (-^ ,2 (l-x)y 2 (l-i/) 3 -cu 2 (l- a ;) 2 y(l-y) 3 (y + x(l-y)) 
-+W (1 - x) y (1 - y f (y + x (1 - y) + y (1 - a;)) 
+ i (i - g 2 ) (1 - x) y (1 - y) 3 (y + x (1 - y) + y (1 - x)) 

+g 2 (l-^) 2 y(l-y) 3 (y + x(l-y))) 



(A2a) 
(A2b) 
(A2c) 

(A2d) 

(A2e) 
(A2f) 
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M 



(8) 

vcs 



xJi; (u(y + x(l- y)) - u'y, ml - t (1 - x) y (1 - y) - q 2 x (1 - x) (1 - yf) 

+a ■ qe'* ■ q (1 - y) [y + x (1 - y) - ^) (6 (1 - x) (1 - y) - (y + x (1 - y))) 

xJ% (u(y + x(l- y)) - u'y, ml - t (1 - x) y (1 - y) - q 2 x (1 - x) (1 - y) 2 ) 

+ (>(l-:r)y(l-y) 3 (y + x(l-y)-£) 

W (1 - x) 2 (1 - yf (y + x (1 - y)) (y + x (1 - y) - 

-uu' (1 - x) (1 - y) 3 (y + x (1 - y) + y (1 - x)) [y + x (1 - y) - 

~ (t - q 2 ) (1 -x)(l- yf (y + x (1 - y) + y (1 - x)) (y + x (1 - y) - 

- q ^l- X ) 2 (l-y)\y + X (l-y))(y + X (l-y)-)^ 

x j£ (cu(y + x(l- y)) - u'y, ml - t (1 - x) y (1 - y) - q 2 x (1 - x) (1 - yf) 
-5(l-j/) 



, (A2g) 



a • £ 



x 4' (vy - J (y + x (1 - y)) , ml - t (1 - x) y (1 - y) - q 2 xy (1 - y 
+ ( - uj 2 y (1 - y) 2 - c/ 2 (1 - x) (1 - y) 2 + x (1 - 1/)) 
(1 - y) 2 (y + x (1 - y) + y (1 - x)) + g 2 y (1 - yf 
-\ (q 2 ~t)(l- yf (y + x(l-y)+y(l- x))) 

xj^ (u;y - (y + x (1 - y)) , ml - t (1 - x) y (1 - y) - q 2 xy (1 - y)) 

+a ■ q'e'* ■ q [ ((1 - x) (1 - y) 2 (1 - 6y) + y (1 - y) (y + x (1 - y))) 
x J'l (uy - u' (y + x(l - y)) ,ml -t(l - x)y(l - y) - q 2 xy (1 - y)) 
+ ( -u 2 (l-x) y 2 (l-yf 



-u ,2 (l-xfy(l-yf(y + x(l-y)) 

+uu' (1 - x) (1 - yf y (y + x (1 - y) + y (1 - x)) 

+g 2 (l-a:)2/ 2 (l-2/) 3 

+ ^(t-g 2 )(l-x)j / (l-2/) 3 (j/ + x(l- 1/ )+j / (l-x))) 

x Jq (u>y- (y + x(l - y)) ,m 2 - t (1 - x) y (1 - y) - q 2 xy (1 - y 
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+a ■ qe ■ q 



(l-y)(l-7y) (y-i) 



(9) 

yes 1 



x J 2 ' (wy-w' (y + x(l - y)) ,m 2 - t (1 - x)y (1 - y) - g 2 xy (1 - y)) 

+ (i - yf (v - 1) 

-J 2 (1 - x) y (1 - y) 2 (y + x (1 - y)) (y - 1) 
+w^'y (1 - y) 2 (y + a; (1 - y) + y (1 - x)) (y - 

+gV(i-y) 2 (y-i) 

~2 (V -t)y(l- yf (y - ~) (y + z (1 - y) + y (1 - ar))) 
x 4' (^2/ - (y + x (1 - y)) , m 2 - t (1 - x) y (1 - y) - g 2 xy (1 - y)) 
(d - 1) J2 (- <V - ml -tx(l- x)) 



F 2 



+ ((a;' - uf x(l-x)-tx(l- xj) J (- (a/ - w) 
Here we have used the definitions 

J 2 (^,m 2 ) 
J 6 (w, m 2 ) 



x, ml 



tx (1 — x 



-Lrf 



, m 2 — u; 2 arccos h (d — 4) 

47r 2 v rriir 



d-l 



d 2 + 6d + 5 
1 

+ d 2 + 6d + 5 
1 



(ml - u; 2 ) J (w, m 2 ) - wA T , 

((d + 5) m 2 - 6cu 2 ) J 2 (cu,m 2 ) 



u; 2 (^ 2 - m 2 ) J (w,m 2 ) 



d 2 + 6d + 5 



umi ( 1 + - ) ) A n 



In Eqs. ( |A3a| ), ( |A3b| ) and ( |A3c| ) we have used the same conventions as f2"T|| , 



2r< I L + 

d-4 



m T , 



16vr 2 



167T 2 /i 

1 1 
1- - 

2 



d-4 



In— + O (d - 4) 
(7.E — 1 — ln47r) 



(A2h) 



(A2i) 



(A3a) 
(A3b) 



(A3c) 



(A4) 



where we introduced the Euler-Mascharoni constant, 7^ = 0.557215, and the scale \x in the 
dimensional regularization scheme we use for the evaluation of the integrals. With J[ and 
J" we denote the first and second partial derivative with respect to m 2 , 

() ' ' ' (A5a) 



J'i (u, ml) 



d{ml 
d 2 



d (ml 



-Ji (w,m 2 ) , 
2 Ji (uJ,ml) . 



(A5b) 
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FIGURES 




FIG. 1. Tree diagrams in virtual Compton scatterin; 
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FIG. 2. Loop diagrams in virtual Compton scattering. 
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